The minimum model for the iron-based superconductors 
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A single band t-U-J-i-Jz model is proposed as the minimum model to describe the supercon- 
ductivity of the newly discovered iron-based superconductors R(Oi- x F x )FeAs and RO\- x FeAs 
(R — La, Ce, Sm, Pr, Nd, Gd). With the mean- field approach, it is found that the pairing occurs in 
the d-wave channel. In the likely parameter region of the real materials, by lowering temperature, 
the system enters firstly the d xy superconducting phase with D^-symmetry and then enters the 
time-reversal-symmetry-broken d xy + id x i_ y i superconducting phase with C^-symmetry. 

PACS numbers: 74.70.-b, 71.18.+y 
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It is well known that the most exciting finding after 
cuprates in the family of superconductors is the discovery 
of R(0^ x F x )FeAs [R = La, Ce, Sm, Pr, Nd, Gd). Soon 
after the announcement of the 26K superconductivity in 
-F-doped LaOFeAs by Kamihara et al [l[ , several groups 
discovered a number of compounds with superconducting 
transition temperatures up to 52K|2j, |3|, U, 15J . Interest- 
ingly, even without F, those compounds may show super- 
conductivity by introducing some oxygen vacancies [fj. It 
is argued that this family of superconductors may share 
some common features with the cuprate superconductors 
and a number of models fH \$k Ifl HolTlll [l2T] have been pro- 
posed to account for the mechanism of its superconduc- 
tivity. Band structure calculations IH 14, 15, [lj| show 
that six d bands of totally ten d bands cross the Fermi 
surface but another group claimed that the parent com- 
pound could be a bad metal 17] with quite low density 
of states at the Fermi surface. A sound experiment com- 
bined with computation show that in the parent com- 
pound there is a phase transition from normal metal to 
stripe-type spin density wave (SDW) around 15OK0. 
This SDW state was subsequently demonstrated by the 
neutron scattering experiments [19|, l2fi. A theoretical 
analysis [13] and ab initial computation [1^] indicate that 
the stripe SDW phase is induced by the frustrated spin- 
exchange interactions and the possible pairing symmetry 
is d xy . These observations somehow hint that the anti- 
ferromagnetic fluctuation may play an important role for 
the superconductivity. 

Although the band calculations indicate that six bands 
(hybridized 2d xy , 2d yz , 2d zx ) cross the Fermi surface, the 
carriers very likely come from those bands lying in the 
a-b plane, i.e., those composed of d xy - and d x 2 _ y 2-ov\Ats. 



However, the numerical result |2l| strongly suggests that 



may qualitatively describe the mechanism of the super- 
conductivity in this family of compounds. Therefore, in 
this Letter, we propose the following Hamiltonian as the 
minimum model to account for the superconductivity in 
electron doped ROFeAs: 

tr,<i,j> j 

+ji s 1 -s 3 + j 2 &-^> w 

<i,j> <<i,j>> 

where a indicates the spin indices; Cj a (Cj !<T ) is the cre- 
ation (annihilation) operator of electrons; rij a is the par- 
ticle number operator and < i,j > and << i,j » indi- 
cate nearest neighbor and next nearest neighbor, respec- 
tively; U describes the Hubbard repulsion and > 
represent the exchange constants. The recent ab initial 
calculation [l6| suggested that J 2 is considerably larger 
than Ji due to the superexchange processes through As 
atoms. 

We note that the Hubbard repulsion U is relatively 
small in the FeAs compounds compared to those in the 
cuprates. The correlation effect is not very strong and the 
FeAs superconductors may share some common feature 
to the so-called gossamer superconductivity [zj l23j |. U 
may have two major effects on the superconductivity: 
One is to expel the s-wave pairing 24] and the other is to 
renormalize the band width and the coupling constants 
Ji^lHl- In such a sense, instead of dealing with model 
(1) we study the following Hamiltonian but keep in mind 
that t and J\ i are the renormalized constants: 



the d x 2_ y 2 band is far above the Fermi surface and d xy - 
orbit indeed contribute a major weight to the Fermi 
pocket. The quasi one-dimensional d yz and d zx bands, where 
though may extend slightly in the iron plane through 
hybridization with p-orbits of As, are unlikely to dom- 
inate the superconductivity. The above arguments di- 
rectly lead to the hypothesis that a single band model 



Heff = t ClC j<a -JiJ2 b hhj ( 2 ) 
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- J 2 E b lAi 
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(3) 



describes the electron pair operator. 
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Mean-field Hamiltonian. Because of positive J\ >2 in the 
system, the electrons favor to form spin-singlet pairs. By 
introducing the <i-wave pairing order parameters d\ =< 



> = 



< h 



>, d 2 =< h 



i,i+x+y 



>= - < 
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tonian: 



> , we obtain the the following mean-field Hamil- 



Hm f = E(k)(ala k + aj k a-k) 

k 

+ J2^ k -E) + 2J 1 d 2 1 N + 2J 2 d 2 2 N, 



(4) 



with N the number of sites and 



E ( k ) = + |A lfc + A 2fe | 2 , 

£ fc = ~2t(cosk x + cosky) - 

A lk = 2Ji<ii7 fc , 

A 2 fe = A 2 = 2J 2 d 2 r] k , 



(5) 



cosky , 77/5 = 2 sin k x sin k y and /i 



where jf. = cosk x 
denotes the chemical potential. The Gibbs free energy at 
a finite temperature T then reads 
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FIG. 1: (color online) The mean-field phase diagram of the 
ground state for n = 1.1. The red lines indicate the phase 
boundaries of d x 2_ y 2, d x 2_ y 2 + id xy and d xy pairing phase. 



G(T) = -2K B T^ln[2cosh(^%)] 



k 



'2K B T' 
& + 2 JidjN + 2J 2 d\N. 



(6) 



For a given density of electron number n, the order pa- 
rameters can be determined self-consistently by the fol- 
lowing equations: 
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(7) 



Ground State. In the ground state, the above self- 
consistent equations are reduced to: 
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FIG. 2: (color online) Finite temperature phase diagram for 
t = 2J and J 2 = 0.9J. 
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Obviously, there are several sorts of solutions for the gap 
functions. If d 2 = 0, we get the conventional d x 2 -y 2 
phase in the ground state; If d\ = 0, we have the d xy 
phase which is isomorphic to the d x 2_ y 2 phase by rotating 



the k space with an angle of 7r/4; For both d\^ 2 =^ 0, there 
is an extra free parameter 6 (exp(i#) = Ji\J 2 \l\Ji\J 2 ) 
which can not be determined by the self-consistent equa- 
tions but can be fixed by the lowest energy. Our nu- 
merical solution shows that only three superconducting 
phases are possible with the variation of t and Ji,2, i.e., 
the d x 2_ y 2 phase for larger J\j J 2 , the d xy phase for quite 
smaller J\j J 2 and the d xy + id x 2_ y 2 phase with 8 = ir/2 
in the intermediate parameter range. The third phase is 
quite interesting because in it the time-reversal symme- 
try is broken. Its point group symmetry is also reduced 
to Cih from D^h of the d xy state. Of course there are 
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also other kinds of mixed solutions such as d xy + d x 2_ y 2 , 
but the corresponding energy is always higher than that 
of d xy + id x 2_ y 2. For large enough t, d\ t 2 — and the 
system is in the normal metallic phase. The mean-field 
phase diagram of the ground state for n = 1.1 and fixed 
J = \J J\ + J i is depicted in Fig.l. In the real mate- 
rials, the band width could be strongly reduced by the 
Hubbard repulsion and the value of J 2 is around 2 J\ or 
even larger as suggested by the ab initial calculation [Hi]. 
In this sense, the superconducting ground state must be 
either the d xy + id x 2_ y 2 or the d xy paired state. 

Finite Temperature Phase Diagram. For given t and 
Ji,2, the thermodynamic phase diagram can be derived 
from the self-consistent equations by taking the order 



parameters tending to zero. Fig. 2 shows the T — n[25| 
phase diagram for t = 2 J and J2 = 0.9 J(Here n does not 
mean the true carrier density in the present mean-field 
approach because we ignore U. In the real materials T c 
around n — 1 could be depressed heavily by U). It is 
found that there are two superconducting phases. One is 
the d xy phase which breaks the U{1) gauge symmetry and 
the other is the d xy + id x 2_ y 2 phase which breaks both 
the U(l) gauge symmetry and the time-reversal symme- 
try. The d xy state has a nodal gap function which allows 
gapless excitations and therefore gains entropy at finite 
temperatures, while the d xy +id x 2_ y 2 state is fully gapped 
with lower energy. The existence of two superconduct- 
ing phases is due to the competing effect between energy 
and entropy. It is emphasized the second thermodynamic 
phase transition could be obtained from the specific heat 
measurement. 

Density of States. One of the important quantities is 
the low temperature density of states, which can be de- 
tected directly by the local probe tunnelling experiments. 
In our d xy + id x 



2 case, it reads 



-y 



2tt e y / U j 2 - |A 1 | 2 cos 2 (2v?) - \A 2 \ 2 sin 2 (2ip) 



,(9) 



where Nf is the density of states of the normal phase 
at the Fermi surface; Ax,2 — 2J1.2d1.2- The numerical 
result of p(uj) is depicted in Fig. 3. Unlike that of the 
d x 2_ y 2 superconductors, there is no node in this time- 
reversal-symmetry-broken superconductor. However, the 
gap function is strongly anisotropic with a small mini- 
mum gap of I Ai| in the spectrum. 

In conclusion, a minimum model to account for the 
mechanism of the iron-based superconductors is pro- 
posed. With the mean-field approach, it is found that 
the most likely superconducting ground state is the d xy 
or d X y + id x 2_ y 2 type. At finite temperature, the d xy 
superconducting state appears first and a second phase 
transition into the d xy + id x 2_ y 2 superconducting state 
occurs in some doping region by lowering the tempera- 
ture. We emphasize that though the real systems are 
multi-band superconductors, we believe that our single 
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FIG. 3: (color online) The left figure is the density of states 
of the d x 2_ y 2 + id xy ground state. The right one is the order 
parameters |di,a| with the variation of J2/J 



band model captures the central physics for the mecha- 
nism of the superconductivity. The multi-band structure 
may only induce multi gaps or renormalizaion of the crit- 
ical temperature without affecting the mechanism. 
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